Focus: A Multifaceted Battle Against Cancer: Extracting Knowledge from Chemical Imaging Data Using Computational Algorithms for Digital Cancer Diagnosis

Fourier transform infrared (FTIR) spectroscopic imaging is an emerging microscopy modality for clinical histopathologic diagnoses as well as for biomedical research. Spectral data recorded in this modality are indicative of the underlying, spatially resolved biochemical composition but need computerized algorithms to digitally recognize and transform this information to a diagnostic tool to identify cancer or other physiologic conditions. Statistical pattern recognition forms the backbone of these recognition protocols and can be used for highly accurate results. Aided by biochemical correlations with normal and diseased states and the power of modern computer-aided pattern recognition, this approach is capable of combating many standing questions of traditional histology-based diagnosis models. For example, a simple diagnostic test can be developed to determine cell types in tissue. As a more advanced application, IR spectral data can be integrated with patient information to predict risk of cancer, providing a potential road to precision medicine and personalized care in cancer treatment. The IR imaging approach can be implemented to complement conventional diagnoses, as the samples remain unperturbed and are not destroyed. Despite high potential and utility of this approach, clinical implementation has not yet been achieved due to practical hurdles like speed of data acquisition and lack of optimized computational procedures for extracting clinically actionable information rapidly. The latter problem has been addressed by developing highly efficient ways to process IR imaging data but remains one that has considerable scope for progress. Here, we summarize the major issues and provide practical considerations in implementing a modified Bayesian classification protocol for digital molecular pathology. We hope to familiarize readers with analysis methods in IR imaging data and enable researchers to develop methods that can lead to the use of this promising technique for digital diagnosis of cancer.     

Isomeric lipid signatures reveal compartmentalized fatty acid metabolism in cancer

The cellular energy and biomass demands of cancer drive a complex dynamic between uptake of extracellular FAs and their de novo synthesis. Given that oxidation of de novo synthesized FAs for energy would result in net-energy loss, there is an implication that FAs from these two sources must have distinct metabolic fates; however, hitherto, all FAs have been considered part of a common pool. To probe potential metabolic partitioning of cellular FAs, cancer cells were supplemented with stable isotope-labeled FAs. Structural analysis of the resulting glycerophospholipids revealed that labeled FAs from uptake were largely incorporated to canonical (sn-) positions on the glycerol backbone. Surprisingly, labeled FA uptake also disrupted canonical isomer patterns of the unlabeled lipidome and induced repartitioning of n-3 and n-6 PUFAs into glycerophospholipid classes. These structural changes support the existence of differences in the metabolic fates of FAs derived from uptake or de novo sources and demonstrate unique signaling and remodeling behaviors usually hidden from conventional lipidomics.     

A four-dimensional computed tomography comparison of healthy and asthmatic human lungs

The purpose of this study was to explore new insights in non-linearity, hysteresis and ventilation heterogeneity of asthmatic human lungs using four-dimensional computed tomography (4D-CT) image data acquired during tidal breathing. Volumetric image data were acquired for 5 non-severe and one severe asthmatic volunteers. Besides 4D-CT image data, function residual capacity and total lung capacity image data during breath-hold were acquired for comparison with dynamic scans. Quantitative results were compared with the previously reported analysis of five healthy human lungs. Using an image registration technique, local variables such as regional ventilation and anisotropic deformation index (ADI) were estimated. Regional ventilation characteristics of non-severe asthmatic subjects were similar to those of healthy subjects, but different from the severe asthmatic subject. Lobar airflow fractions were also well correlated between static and dynamic scans (R² > 0.84). However, local ventilation heterogeneity significantly increased during tidal breathing in both healthy and asthmatic subjects relative to that of breath-hold perhaps because of airway resistance present only in dynamic breathing. ADI was used to quantify non-linearity and hysteresis of lung motion during tidal breathing. Non-linearity was greater on inhalation than exhalation among all subjects. However, exhalation non-linearity among asthmatic subjects was greater than healthy subjects and the difference diminished during inhalation. An increase of non-linearity during exhalation in asthmatic subjects accounted for lower hysteresis relative to that of healthy ones. Thus, assessment of non-linearity differences between healthy and asthmatic lungs during exhalation may provide quantitative metrics for subject identification and outcome assessment of new interventions.     

Evaluation of four-dimensional cone beam computed tomography ventilation images acquired with two different linear accelerators at various gantry speeds using a deformable lung phantom

We evaluated four-dimensional cone beam computed tomography (4D-CBCT) ventilation images (VI_CBCT) acquired with two different linear accelerator systems at various gantry speeds using a deformable lung phantom.The 4D-CT and 4D-CBCT scans were performed using a computed tomography (CT) scanner, an X-ray volume imaging system (Elekta XVI) mounted in Versa HD, and an On-Board Imager (OBI) system mounted in TrueBeam. Intensity-based deformable image registration (DIR) was performed between peak-exhale and peak-inhale images. VI_CBCT- and 4D-CT-based ventilation images (VI_CT) were derived by DIR using two metrics: one based on the Jacobian determinant and one on changes in the Hounsfield unit (HU). Three different DIR regularization values (λ) were used for VI_CBCT. Correlations between the VI_CBCT and VI_CT values were evaluated using voxel-wise Spearman’s rank correlation coefficient (r).In case of both metrics, the Jacobian-based VI_CBCT with a gantry speed of 0.6 deg/sec in Versa HD showed the highest correlation for all the gantry speeds (e.g., λ = 0.05 and r = 0.68). Thus, the r value of the Jacobian-based VI_CBCT was greater or equal to that of the HU-based VI_CBCT. In addition, the ventilation accuracy of VI_CBCT increased at low gantry speeds.Thus, the image quality of VI_CBCT was affected by the change in gantry speed in both the imaging systems. Additionally, DIR regularization considerably influenced VI_CBCT in both the imaging systems. Our results have the potential to assist in designing CBCT protocols, incorporating VI_CBCT imaging into the functional avoidance planning process.     

In vitro dynamics of chromatin organization and migration

The organization of eukaryotic chromatin is not static but changes as a function of cell status during processes such as proliferation, differentiation, and migration. DNA quantification has not been used extensively to investigate chromatin dynamics in combination with cellular migration. In this context, an optimized DNA-specific, nonperturbant method has been developed for studying chromatin organization, using the fluorescent vital bisbenzimidazole probe Hoechst 33342: this property has been described by Hamori et al. (1980). Computer-assisted image analysis was used to follow migratory activity and chromatin organization of L929 fibroblasts during in vitro wound healing. Cell movements were analyzed using an optical flow technique, which consists in the calculation of the velocity field of cells and nuclear movements in the frame. This system allows the correlation of cell migration and position in the cell cycle. It makes it possible to study chromatin dynamics using a quantitative analysis of nuclear differentiation reorganization (nuclear texture) and to correlate this with migration characteristics. The present system would be of interest for studying cell-extracellular matrix interactions using differing substrates, and also the migratory response to chemotactic factors. Such a model is a prerequisite for gaining better understanding of drug action.     

Dream Imagery and Emotion.

The relationship between prominent visual imagery and emotion within dreams was investigated in relation to E. Hartmann's (1996) contextualizing image (CI) theory and M. Seligman and A. Yellen's (1987) dual imagery theory. Fifty-nine students recorded dreams over a 2-week period and submitted 115 dreams for analysis. Participants recorded ratings of emotion type and emotion intensity in each scene. Prominent visual images were identified and scored for intensity and detail by independent judges. As hypothesized from Hartmann's theory, there was a significant positive relationship between CI intensity and emotion intensity in the CI scene, emotion intensity generally peaked in the CI scene, and dreams containing a CI had higher overall ratings of emotion intensity than non-CI dreams. The result for the correlation of detail of prominent imagery with emotion was inconclusive, with a low positive correlation across CI scenes. This raises the possibility that the CI is not a unitary construct. (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

The biometry of gills of 0-group European flounder

The gill surface area of 0-group, post-metamorphic Pleuronectes flesus L. was examined using digital image analysis software and expressed in relation to body mass according to the equation log Y=loga+c logW ( a =239·02; c =0·723). The components that constitute gill area, total filament length, interlamellar space and unilateral lamellar area were measured. The measurement of the length of every filament on all eight arches showed that commonly used methods of calculation can lead to an under-estimation of up to 24% of total filament length. Direct measurements of unilateral lamellar area with digital image analysis showed that previously reported gill area data for the same species was over-estimated by as much as 58%. In addition, in this species the neglect of gill pouch asymmetry after metamorphosis, can bring about a 14% over-estimation of total gill area.     

Non-invasive image analysis evaluation of growth during plant micropropagation

Microcomputerized video image analysis was adapted for rapid, objective, and non-intrusive quantification of shoot growth and development for plants growing in vitro. Custom-developed staging arrangements were essential to insure accurate viewing and representation of the plants in each of three standard culture vessels. Shoot length measurements from digitized culture images were strongly correlated with length measured manually ex vitro. Image analysis weighted density measurements of proliferating microcultures (even with irregular growth habits) provided a reliable indicator of shoot culture fresh weight. Non-destructive time course evaluations of growth rate and quality were demonstrated.Abbreviations FW fresh weight - WD weighted density     

Total root length as estimated from small sub-samples

A method is proposed for estimating the total length of a root system from sub-samples. This method is based on the measurement of the length and diameter of small pieces of roots, and on the measurements of the bulk density of root sub-samples. It is assumed that roots are cylinders with a given bulk density. The length and diameter of small root pieces are measured by image analysis. A weighted quadratic mean (W.Q.M.) root diameter is then calculated and used in estimating the root length. This W.Q.M. diameter is defined as the real mean diameter of an equivalent single root with the same length and volume as the tested root system. The accuracy of prediction is demonstrated for one theoretical root system. The standard deviation of estimation can be calculated using sampling simulations.     

An inflatable minirhizotron system for root observations with improved soil/tube contact

Commonly used minirhizotrons consisting of a transparent tube inserted into the soil seldom attain good contact between the tube and the soil, which leads to root growth occurring in a gap rather than in the soil. A new system is described involving an inflatable flexible rubber wall, made from a modified motorcycle tube. Pressure ensures a proper tube/soil contact so that the environmental circumstances for root growth along the tube more closely correspond to those in the undisturbed soil. Before the endoscope slide is introduced into the minirhizotron for taking pictures, the inflatable tube is removed, so that there is no-often opaque-wall between the endoscope and the roots. This improves the picture quality and facilitates the analysis of root images.